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Abstract 

Equations of motions and energy-momentum density tensors are 
obtained for a dispersive and dissipative medium sustaining electric 
and magnetic polarizations, using Lagrangian formalisms. A previous 
work on the subject by the authors has been simplified by reduction 
in the Lagrangian of the number of independent vector fields, with 
which the sink modes (associated with the dissipation) interact. The 
paper also formulates the dynamics of the electromagnetic field in the 
medium in a covariant (relativistic) manner. We discuss (and compare) 
the results, especially from the point of view of the positivity of the 
energy density. 
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1 Introduction 

Lagrangian methods provide a systematic and unified approach to dynamic 
phenomena. Every mature subject in physics can and should be described 
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by a Lagrangian. To name a few: particle physics is described by the "stan- 
dard Model" Lagrangian, General relativity by the "Curvature Scalar" La- 
grangian and so on. 

Nevertheless, there are some phenomenological treatments that are not 
obviously expressible by a Lagrangian formalism and have led to ambigui- 
ties and controversies. An instance of this is the dynamics of polarization 
in dispersive-dissipative materials. While the equations of motion (shown 
below in equation (10) and equation (11)) have been in frequent use, the 
derivation of these equations from a Lagrangian has been wanting for some 
time. As a consequence, basic quantities like energy and momentum densi- 
ties, that are integral parts of the Lagrangian method, have become variously 
formulated and subject to dispute. A brief account of the history of energy 
densities is found in [1] . Following that account and adopting a Lagrangian 
formulation, that work has unambiguously derived an expression for the 
energy density. 

Relativistic treatments of dissipative electromagnetic media, but not 
based on a Lagrangian formalism, were given in [2]. A comparison of three 
different approaches for relativistic treatments of dissipative electromagnetic 
media was undertaken in a book by Fano, Chu and Adler [3]. In this book 
the authors compared the approach they endorse with two other approaches, 
denoted the Minkowski and Amperian approaches. The Minkowskian ap- 
proach was criticized by the authors because it assumes a linear scalar rela- 
tion between the displacement field D and the electric field E, and between 
the magnetic field H and the magnetic flux density field B. The authors 
argue that this relation is not invariant under a general Lorentz transfor- 
mation and thus is not objective but subjective and can only be regarded 
as true in a a particular frame of reference. Furthermore, even according 
to an observer in rest with respect the electromagnetic medium, the lin- 
ear scalar relations are not generally true since some magnetic materials 
have a hysteresis curve and thus B and H are not linearly related, while 
other materials possess a tensorial dielectric constant instead of a scalar 
one. Fano, Chu & Adler also claim that the Lorentz force resulting from 
the Minkowskian formalism contradicts experimental evidence regarding po- 
larized and magnetized materials. The Amperian approach is criticized in 
the book on the basis that the resulting Lorentz force contradicts experi- 
mental evidence regarding magnetized materials. The authors of this paper 
do not follow the approach of Fano, Chu & Adler, since it is not based on 
a variational approach and thus one cannot canonically derive an energy 
expression. Furthermore, Fano, Chu & Adler do not regard the magneti- 



2 



zation and polarization of a material as degrees of freedom subject to their 
own dynamics; rather, these are treated as given source terms for the the 
electric and magnetic fields. It should be noted however, that both the en- 
ergy density obtained by them and that using the Amperian approach are 
positive. 

Other approaches, notably for physical systems containing both vis- 
cous flows and electromagnetic fields are discussed by Felderhof [4] and 
Rosensweig [5]. The first work postulates an energy-momentum tensor 
(equation (3.12) in [4]); in our formalism we will aim at deriving the en- 
ergy momentum tensor from a Lagrangian using the well known canonical 
formalism (see below). The second work [5] is based on an "engineering" 
type of approach. Rosensweig has postulated balance equations in a moving 
frame for linear momentum, angular momentum and energy and has derived 
"constitutive equations" using the first and second laws of thermodynam- 
ics. The energy momentum tensor was not considered in this work, nor a 
Lagrangian to derive the energy momentum tensor from. 

It is evident that the energy density, written as T$, is a component of 
the energy momentum tensor Tj ([9]- [12]). This is uniquely derivable from 
the Lagrangian density £(«&, Ukj), which is composed of the field variables 
Uk and of their derivatives Ukj- (A symbol after the comma represents a 
derivative.) The formal definition of the energy- momentum density is given 
by 

Summation is implied for repeated indexes. In our convention the indexes 
take the value for the time component and the numbers (1,2,3) for the 
remaining, space components. The dot over a quantity will signify a time 
derivative. 

The expression for the energy found in [1] was positive. The essential idea 
there was the introduction of added degrees of freedom (named sinks), which 
are depositories of the the energy residing in the electromagnetic field and in 
the polarization of the medium. Physically, the sinks are macroscopic modes 
of motion in the solid (and less obviously, in its surrounding) different from 
the polarization modes, but weakly interacting with them. These degrees of 
freedom have been given microscopic interpretation in the thermodynamic 
limit, in terms of a set of interacting oscillators with a special interaction 
matrix [6]. The need to include additional degrees of freedom to treat the 
irreversible behavior of polarizable systems was noted in a classic text [7]. 
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The proof that any set of differential equations can be embedded in a larger 
set that can be derived from a variational principle was supplied by Bateman 
[8]- 

Our aim in the present work is twofold: first we simplify our previous 
formulation by reducing the number of vector degrees of freedom. Thus, the 
formalism presented here contains two vector degrees of freedom less than 
the previous one. The resulting Lagrangian is considerably simpler without 
compromising the power of our formalism. Secondly, in later sections, we 
give a relativistic covariant formalism, making the theory valid also in a 
moving system of reference. The motivation for this is twofold: 

1. It is well known that any physical theory should obey the laws of 
special relativity and should be invariant under Lorentz transformation 
- in this paper we construct a theory along the same lines as our 
previous theory [1], which satisfied these requirements. 

2. Our previous theory did not contain spatial derivatives of certain quan- 
tities; the principle of covariance has helped us put those derivatives 
in. 

Further, we discuss the implications. Some of the equations in the earlier 
sections have already appeared in [1]; however, they are indispensable both 
for a self-contained reading of the simplified theory and as an introduction 
to and comparison with the covariant theory. 



2 Time independent polarizations 

Our first task is to formulate the Lagrangian density of the electromagnetic 
and polarization and magnetization fields for a material whose polarization 
are static. We use CGS Gaussian units. Connections between the electric 
displacement field D and the electric field E and polarization P, and, like- 
wise, between the magnetic field H and the magnetic induction field B and 
magnetization M are then provided by the relations: 

D = E + 4ttP. (2) 

H = B — AttM. (3) 
The above fields satisfy both the homogeneous Maxwell's equations: 

curlE + -B = 0, divB = (4) 
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and the inhomogeneous equations: 



1 u 4tt - 

curlH D = — J, divD = Ann (5) 

c c 

We further make the assumption that both free charges p and currents J 
are zero. 

It is not possible to derive these equations from a Lagrangian density 
expressed in terms of the fields. However, this difficulty is traditionally 
overcome by representing the fields in terms of vector A and scalar <J> po- 
tentials, as follows: 

E = -V$ - ^A, B = curlA (6) 

Using these definitions one can see that the homogenous equations (4) are 
satisfied automatically. The inhomogeneous equations (5) are the outcome 
of the functional derivation of the Lagrangian below: 

£ = &EB + ^pe + £mh 
C EB = ^[E 2 -&] 

C PE = PE 

Cmh = M-H (7) 

In forming the derivatives it is understood that all field quantities as given in 
terms of the potentials of equation (6). The energy density is then obtained 
from equation (1): 

T° = ?t.A-£. (8) 
dl 

After some algebra this yields for the Lagrangian of equation (7) the result: 

T ° = — [E 2 + H 2 + (4ttM) 2 ] (9) 

The last term M 2 can be usually discarded for static magnetization, since it 
represents nothing but a time-independent constant. However, we include 
it for reference in future sections, where the electric polarization and the 
magnetization will become dynamic degrees of freedom. 
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3 Equations of motion with dissipative terms 



The equations that treat the development of the polarization P by an applied 
electric field E, and the magnetization M as induced by an external magnetic 
field H are, as appearing in e.g., [13], the following: 

P + T e P + u?P = u$E (10) 

M + T h M + loqM = FluqH (11) 

with all symbols and notation as defined in [13]. These equations hold for 
t > 0. Fourier transformed to the frequency (a;) domain, the equations for 
the frequency components of the fields are written, after a transposition of 
terms as 

[-uj 2 - iuT e + uj 2 r }P{uj) = u 2 p E{uj) + P(t = 0) - {icu - T e )P(t = 0) (12) 

and similarly for the magnetization. (Here time-derivatives of the polariza- 
tion were integrated by parts in the manner: / °° P(t)e tujt dt = [P(t)e tujt ]^ — 
iuj Jq°° P(t)e luJt dt and neglecting the field at infinity, etc. As discussed in [7], 
there can be a polarization at t = due to permanent dipoles.) One can 
introduce the electric and magnetic susceptibilities, defined by 

oSl/u 2 

XcM = V-JL ^ (13) 

Xm{Uj) = 1 _(Jp_ i *T h ( 14 ) 
to write the equations of motion for the frequency components as: 

P(o)-(iLu-r e )P(o). 



UJ 2 



M(w) = Xm ( w )[ff( w ) + — U v v v ; ] (15) 

The electric and magnetic susceptibilities are complex functions with related 
imaginary (3 s ) and real (3ft) parts. When the electric losses are sufficiently 
low T e << 1, the electric susceptibility can be approximated in the vicinity 
of the resonance by: 

3[Xe] = 1 K[ Xe ] = -^^[Xe] (16) 

2 
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Analogous relations can be derived for the real and imaginary parts of the 
magnetic susceptibility for « 1: 

^Xm] = ^ J^ W — , nXm}=-^Z[ Xe ] (17) 

2 

Returning now to equation (10) and equation (11), for self-consistency these 
have to be supplemented by the Maxwell-equations given in the previous 
section and derived from the Lagrangian of equation (7) that contain the 
contravariant four- vectors A a = (A, $>). 

To take account of the dissipative nature of the processes we introduce 
a set of new fields, designated canonical fields. These will be formally dis- 
tinguished from the previous, physical fields by writing them in lower case 
symbols. The introduction of these fields follows theories by [14] -[16]. The 
canonical fields (taken to be real) are made to be part of a Lagrangian for- 
mulation. We note that for the preceding physical quantities (written in 
upper case symbols), this cannot be achieved [14]. 

To make our new notation explicit, we shall work with electric polariza- 
tion p a and magnetization m a and the electromagnetic fields e a , h a and the 
(contravariant) vector-potentials a a . It will be later demonstrated that these 
new, canonical fields contain a time dependent factor which (partly, but not 
fully) neutralizes the dissipative process. Otherwise, the choice of the fields 
(and of the various constant factors) is guided by the requirement that we 
recapture the form of the equations of motion and the energy density cur- 
rently widely employed in solid state optics in the appropriate limiting case 
of no dissipation. (This requirement impairs at times the formal symmetry 
between the modes linked with the polarization and the magnetization, re- 
spectively.) The physical meaning of the canonical fields will be made clear 
later by comparing their equations of motion with those of the physical 
fields. (For a similar procedure see [14].) 

In addition, we shall introduce two new scalar fields, the so called "sink" 
fields: the first being VP, that incorporates macroscopic degrees of freedom 
associated with the dissipation of the polarization (analogous to the degree of 
freedom denoted with the same symbol in [16]) and the second, the sink- field 
<I> which accounts for the dissipation in the magnetic mode (and originates 
physically in atomic spin or angular momenta in the solid). 
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4 The Lagrangian with sink terms 



The addition of sink fields requires the extension of the Lagrangian as fol- 
lows: 

£ = £-eb + C-pe + £>mh + C p q, + + £<t> + £<!> (18) 

In this Lagrangian we have first the Lagrangian of equation (7), there writ- 
ten in terms of the electromagnetic fields, and now rewritten in the new 
lower-case variables as C e b + C pe + C m h- The above expression also con- 
tains the Lagrangians in the polarization modes, the magnetization modes, 
the electromagnetic fields, the ^ and <I> sink fields and then the interaction- 
Lagrangians between polarization and electric field, magnetization and mag- 
netic field and, finally, the polarization and magnetization with their respec- 



tive sink fields. The various terms are written out as: 

C pe = pe! (20) 

C m h = mhi (21) 

£ p * = ^|(P 2 -K 2 -*V-2PP*) (22) 

£ m $ = -^(m 2 -(w 2 -$ 2 )m 2 -2mm$) (23) 

It U!q 

U = ^etf 2 -^* 2 (24) 

U = Ia6 2 -^$ 2 (25) 



having chosen the x (or 1) the direction of the polarization and mag- 

netization vectors . In the above, we have expressed and as harmonic 
modes of motion (when there is no coupling to other fields.) This choice 
seems to allow the simplest type of representation of unspecified (and, so 
far, arbitrary) degrees of freedom. The symbols k, p, A, v are constant, non- 
negative coefficients, whose values depend on the nature of the relaxation 
mechanisms. 

We next derive C with respect to each degree of freedom. In accordance 
with the Euler-Lagrange equations, this procedure leads to the equations of 
motion. In detail, for the variables p and m the equations of motion take 
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the following form: 

p - (* + ^ 2 - uj 2 r )p = u? v e x (26) 
m - ($ + $ 2 - wg)m = Fu^/ii (27) 

The electromagnetic equations of motion are identical to Maxwell's equa- 
tions in the new (small letter) variables. Lastly, for the sink coordinates one 
finds 

$ + * 8 fr0-P*>) (28) 



re kuj. 



dt 



■■ v , 1 d(m(m — m$)) , . 

* + a* = a^ « (29) 

The expression for the energy density Tq is found from equation (1), 
namely, 

^EJI-W (30) 

with the summation over all degrees of freedom. We separate the electro- 
magnetic, the polarization and the magnetization parts 

T 0° = ( T o)e,h + (To)p^ + (To) m ,<l> (31) 

where the first term has the well known form of the electromagnetic energy 
density, already given in equation (9): 

(T °) e ,/ l = ^[e 2 + ^ + (47rm) 2 ] (32) 
The others take the forms: 

(T °k* = ^[(p-p^+^ + f^ + f^ 2 (33) 

1 ^i! (34) 



2u 2 p dt 



and 



(T °° )m '* = ^2[(™-^) 2 +^ 2 ] + ^ 2 + ^ 2 (35) 
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1 I • 2 i 2 2\ i ^/k2 , V ^2 , ^ 2A2 

1 dm 2 ■ . . 

- 2Kjir <t < 36 > 

The first- written forms are sums of squares with non- negative coefficients; 
this ensures that each part of the energy density is positive (non-negative). 
The usual expressions for the energy density, e.g. in [12] or [13], differ from 
the above by the presence of the terms in $ and VP and their time derivatives. 
In the following section we shall eliminate these variables by making use of 
the equations of motion, equation (28) and equation (29). 



5 A simple solution 

To work out a fully solvable case, we postulate that the coefficients p and v in 
equation (24) and equation (25) vanish. Then the equations of motion for the 
sink-modes, equation (28) and equation (29), can be integrated. We shall 
carry through the electric polarization case, but an analogous development 
holds for the magnetization sink- variable. We assume the following initial 
conditions for \P(t) 

(0) = 0, (0) = T e /2 (37) 

With these choices the differential equation (10) is regained for short times, 
t << y~, as will be shown presently. (This is similar to the procedure in 
[14]). Then, from equation (28), 

9- 1 9 \ 2 (t) 1 d{p2{t) * ] (38) 

*-^ p w p{t) ~^4 dt (38) 



A first integration gives 



which leads to 



2nuip at nujp 



klj£ 2kuj* ot 



where the constant is 



c^a+^-mm (41) 

2 KU£ KUJ* 
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Integrating once more and fitting the constant so as to satisfy the first initial 
condition in equation (37), we finally obtain: 

1 V + 



p 



2 W^f) J 
+[ r, p(Q)p(Q) 1 /* g (42) 



P •> U ± "1 



One can now substitute this expression and equation (39) into equation 
(34) to obtain, after considerable simplification, an expression for the energy 
density arising from the time varying "canonical" polarization pit) in the 
form: 

n 1 /P 2 (t) + klu 2 C 2 „ o \ 

( T o°W = ^|( ~ Y^gf + «rV(*)) (43) 
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In the case of slow relaxation, when T e is small (quantitatively, when T e ^fE « 
1), then so is p(0). (This will be confirmed in the next section.) This results 
in the quantity C 2 being a small quantity of the second order, which can be 
neglected: 

^ = M^m + ^ p2{t) ) (44) 



For the part of the energy density involving the " canonical" magnetization 
m(t), the variables n(t) and $(t) having been eliminated through their equa- 
tions of motion, a similar procedure gives 

(ToV, = ^(-5Sw + ^ 2 W) (45) 

The above two expressions for the energy densities, equation (44) and equa- 
tion (45), resemble those in equation (11) of Ruppin [13], except that they 
are written in the canonical (small letter) variables, rather than in the phys- 
ical variables (the relations between these will be presently obtained), and 
that in the denominators they contain the polarization fields squared. The 
role of these is to damp out "kinetic" energy in the polarization motion, 
associated with the motion of charged (or spinning) matter. 

The main results of this section, equation (43) and its analogue for the 
magnetization energy density, are exact and contain non-perturbative cor- 
rections to the energy density, due to the presence of the sink degrees of 
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freedom. While exact, they are model dependent in the sense that sinks rep- 
resented by different Lagrangians would lead to different energy densities. 
This is clear, due to the presence in the energy densities of the parameters k 
and A that were introduced in the Lagrangian in equation (24) and equation 
(25). This outcome was anticipated some time ago in [17]. (We also note 
the opposing view in [18].) It is of interest to note that the non-dissipative 
limit is not regained when T e ,Th — > 0, but only when also n,\^oo. It 
can be shown that in these limits the canonical fields are identical to the 
physical fields. 



6 Interpretation of the fields 

We now find the relation of the canonical fields to the physical fields in the 
presence of dissipation and sink modes. We regain the original equations 
of motion, equation (10) and equation (11), for the physical polarization 
variable, as follows: We postulate 

p(t) = e*C'P(t) (46) 
m(t) = e* W M(t) (47) 

This turns (the vector form of) equation (26) into the following: 

P + 2V(t)P(t) +w r 2 P = ^e-*( t 'e (48) 

and likewise for the magnetization variables. Then, from equation (37), for 
short times < t << 

tf(i) « *(0) = T e /2 (49) 

We then obtain 

P + T e P(t)+uj 2 r P = LU 2 p e-^e (50) 

Recalling equation (10) and equation (11), we can thus extrapolate to later 
times so as to identify 

e(t) = e*®E{t) (51) 
h{t) = e*®H(t) (52) 

and regain equation (10) and equation (11). This provides a physical mean- 
ing for all the "canonical" variables, as those fields in which the decay has 
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been reinstated. On the other hand, the decay is itself dependent on the 
fields. (Cf. [19].) Furthermore, the Maxwell equations for the physical fields 
are also modified in the dissipative-polarizable medium. This feature (of a 
modified Maxwell equation) also appears in [15] and [16] (eq. (14) and eq. 
(12), respectively). 



7 A covariant dissipative Lagrangian density 

In this section we introduce a covariant Lagrangian density which is a 
Lorentz invariant generalization of the Lagrangian density described in equa- 
tion (18). Again we introduce a set of new fields, which are formally dis- 
tinguished from the previous, physical fields by writing them in lower case 
symbols. Explicitly, we shall work with the scaled polarization and elec- 
tromagnetic tensors p a @ and f a @ and other lower case quantities that were 
introduced earlier. Notice that the Lagrangian introduced here will contain 
less terms since both the magnetization and the polarization are now part 
of the polarization tensor p a/3 . This Lagrangian has the following parts: 

C = Cf p + C p \& + (53) 

In which Cf p is defined in equation (89) in the appendix (although an up- 
per case symbols to lower case symbols transformation is needed). £ p y is 
obtained by generalizing equation (22) as follows: 



d&^&Vofi ~ H ~ d^d^) PafS p al3 - 2p af ,d^d^] (54) 



Cp^, can be written as a sum of two terms: one depending on the magneti- 
zation, while the other depends on the polarization. 



4* = -i[-4(^p) 2 + (^p) 2 + K-^(W 2 + (^*)V 

ujZ c c 
2 

+ — p-pV -2p-dipdiV] 



C r m = -^\--{d t mY + (d i mY + (^.--(d t ^y + (v^Y)m 



2 • • 

+ — m • m^> — 2m ■ dirndl] (55) 
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Finally is obtained by generalizing equation (24) as follows: 

= \[*d^d^ - p* 2 ] = I[£(l(^) 2 - (W) 2 ) - p^ 2 } (56) 

The field equations with respect to p a/3 , ^ and the electromagnetic fields can 
be obtained by taking the variational derivative of C given in equation (53) 
with respect to p a/3 ,^ and a a . The electromagnetic field equations were 
already derived by taking the variational derivative of equation (89) with 
respect to a a and further change is unnecessary (except for that from upper 
to lower case symbols ). Taking the variational derivative of C with respect 
to ^ leads to the equation: 

kd^ + p* = \d^{d^p a p - d,f!lp a p)] (57) 

Taking the variational derivative of C with respect to p a p leads to the equa- 
tion: 

d^p^ + K 2 - d^d^ - d^) Pa/3 = -\ujjf a p (58) 

In the "homogeneous" case in which the spatial derivatives of both the 
polarization tensor p alS and the scalar field * vanish, the tensor equation 
(58) can be written as two vector equations: 

P- (* + ^ 2 - u 2 .)p = u 2 p e (59) 
rh - (* + # 2 - u? r )m = -u? p b (60) 

In which lo 2 . = c 2 uj 2 and lo 2 = ^c 2 lo 2 . We can clearly see that equation (59) 
is the same as equation (26). As for equation (60) by using the equality (3) 
this can be rewritten as: 

fh - + 4> 2 - u)l)m = -uj 2 p h (61) 

In which uj 2 , = u) 2 + 47ro; 2 . This equation can be identified with equation 

•J 1 

(27) by choosing F = — | and equating the magnetic and electric dissipative 

modes ($ = *). Notice that in magnetized materials the magnetic field is 
in a direction opposite to the magnetization as can be seen from figure 1, 
which explains the negative sign in the above equation. 
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Figure 1: Magnetic field lines in a box of permanent uniform magnetization 
- a cross section view. The direction of magnetization can be inferred from 
the way the magnetic field lines look outside the box magnet. The normal 
component of the field changes sign across the material boundary, unlike in 
the dielectric case. 

7.1 The energy momentum tensor of the dissipative 
Lagrangian density 

Using the dissipative Lagrangian density defined in equation (53) we can 
calculate the energy momentum tensor by using the formula given in (1). 
Doing so we obtain the following expression: 

= T%[fp] + TfaV] + T^} (62) 

T@[fp] is calculated in Appendix (A. 3) in equation (93). It remains to obtain 
expressions for Tj^[p^] and T^[^]; these are given below: 

+ dPvdaVpTep^-SPCpv (63) 

and 

r£[*] = «^*a a * - s&u (64) 

Let us look at the energy density Tq. For Tq[\I/] we obtain the positive 
quantity 

r [f] = ^[^ 2 + /c(W) 2 + p* 2 ] (65) 
For T$\pV] we obtain 

?o [p*] = ^[(d^m - d^mf - (d^p- d^p) 2 + Q? T {m 2 - f)] (66) 
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This expression contains positive magnetization contributions and negative 
polarization contributions. The positive magnetization contributions can 
overcome the negative magnetization contribution of Tq[/p] given in equa- 

" 2 

tion (94) of Appendix (A. 3), provided that % > 2ir. However, the negative 

polarization contribution to Tq results in the uncertainty whether the total 
polarisation contribution is negative or positive, since this depends on the 
solution of equations (88,57,58). Those equations are difficult to solve since 
they are complicated non-linear partial differential equations. In the follow- 
ing subsection we will present a solution for a particular, simple situation. 
Notice that this situation is completely different from the situation in the 
non-covariant model for which the positiveness of the energy density was 
assured. To conclude this subsection we calculate the contributions to Tq: 

n = T [fp\+T^\+T^\ (67) 

since Tg[/p] given in equation (95) is the Poynting vector, we can consider 
Tq[p^] and Tq[^] as corrections to the Poynting vector. Those are given by: 

T >*] = — [Vp k p k - Vm k m k - p k (V^p k + Wp fe ) 

CLOp 

+ m fc (Wm fc + Wm fc )- V^^{m 2 - f)} (68) 

and 

2*[tf] = — (69) 

One should notice that both these corrections vanish for the homogeneous 
case in which the spatial derivatives vanish. 



7.2 A covariant solvable case 

In this final subsection we would like to elucidate the nature of the energy 
density component Tq by introducing a simple solution of equation (57). 
We assume a homogeneous situation in which all spatial derivative vanish. 
Furthermore, we assume the p = and introduce the constant k = In 
this case we obtain the equation 

1 



= ^d t [m-m-^m 2 -p-p + ^f'} (70) 

From this equation we can obtain an expression for \E', as follows: 

,1 d t (rh 2 — p 2 ^ 
v 2 Ku>^ + m 2 — p 2 
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(In which we have neglected a constant that can be considered small). Using 
equation (71) we can rewrite equation (65) as 




2 " 8" ^/to; 2 + m 2 
To this we need to add the homogeneous energy density 

= i bi [{fh ~ *™ )2 _ ^ 2 + w2( ™ 2 ~ ?2)] (73) 

We will consider to separate cases: 

1. p = 0, m = rax. In this case we obtain 

1 Til 2 

T °[*] + T >*] = — [ w 2 m 2 + (74) 

which is very similar to equation (45). This expression is positive 
and with the correct choice of parameters can overcome the negative 
magnetization contribution of Tq [fp] 

2. p = px, m = 0. In this case we obtain 

-i -2 2 

1 , o;^ 

2 p 2 — KOJ 2 Wp' 

which is very different from equation (44). This expression is positive 
only if p 2 > kuj 2 and if p changes fast enough to overcome the negative 

part — %p 2 ; otherwise, the energy density is negative. 



T^]+^m = -J^— 2 -^P 2 } (75) 



8 Conclusion 

In this paper the Lagrangian method for a dissipative medium (capable 
also of sustaining electric and/or magnetic polarizations) was carried out 
to obtain equations of motion, energy-momentum densities, etc. To apply 
the method it has been found necessary to introduce additional degrees of 
freedom (" sink- modes" ) , associated with decay mechanisms in the electrical 
and magnetic modes. 

Two versions have been formulated: The first one, which is not invariant 
under Lorentz transformation, has led to equations for both the electromag- 
netic fields (a slightly corrected set of Maxwell equations) , polarization and 
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magnetization and, using a model for the sink variables, has unambiguously 
yielded energy densities, shown in equation (43) - equation (45). These are 
positive and resemble the corresponding results in [13]. The results obtained 
are similar to the ones obtained in [1], however, they have been obtained 
here with a much simpler Lagrangian, containing fewer degrees of freedom. 
To be precise, two vector degrees of freedom (or six independent quantities) 
have been omitted from the beginning, thus making the formalism more 
useful for numerical calculations. 

Secondly, a covariant Lagrangian was formulated, one that is invariant 
under the Lorentz transformation and contains spatial as well as tempo- 
ral derivatives. From this, we have derived equations for both the electro- 
magnetic fields, polarization and magnetization . Calculating the energy 
momentum tensor we have derived expression for both the energy density 
and the generalized Poynting vector. However, we find that with a specific 
model for the sink modes, the energy density derived from this covariant 
Lagrangian is not positive definite. This points to the need to introduce a 
better model for the sink degrees of freedom, perhaps with better physical 
underpinning. 

APPENDIX 

A The covariant formalism 

In a covariant formalism of the electromagnetic theory one uses a four dimen- 
sional space-time formalism instead of the usual distinction that is practised 
between spatial and temporal coordinates. In this formalism the invariance 
of the Lagrangian of the electromagnetic theory under Lorentz transforma- 
tion is clear. 

A.l Four dimensional notations 

The four dimensional coordinate are defined as: 

x° = ct, x 1 = x, x 2 = y, x 3 = z (76) 

In this formalism there is a difference between coordinates with upper and 
lower indexes. For example: 

x = ct, x\ = -x, x 2 = —y, x 3 = —z (77) 
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The connection between upper and lower coordinates is established through 
the metric matrix g a p, which is defined as 

Q/ 3 = Q/3 = diag(l, -1,-1,-1) (78) 

(we shall use Greek letters for 4 coordinates according to the well established 
custom) such that 

W 7/3 = ^ = diag(l, 1,1,1) (79) 

Using g a p we can write x a = g a px^ '. This notation will be adapted for any 
four vector 



A — (A , A), A a — {A , A), A a 



da/3 



A? 



For the four dimensional derivatives we will use the notation 



d a = 



,d_ 

dx ' 



-v), d a 



d 

dx° 



(80) 
(81) 



A. 2 Electromagnetic theory in four dimensions 

The vector and scalar electromagnetic potentials are replaced by a single 
four-vector 

A a = ($, A) (82) 

Using equations (6) and (82) we can define the contravariant and covariant 
antisymmetirc tensors: 



F af3 = Qa A f3 _ q/3 





( 


-E x 


-Ey 


-E z \ 




E x 





-B z 


B'y 




EV 


B z 





-B x 




\ EZ 


-By 


B x 


o / 


( 





E x 


EV 


E z \ 




-E x 





—B z 


By 




-E y 


B z 





-B x 


v - 


-E z 


-By 




o J 



(83) 



(84) 



Analog tensors can also be defined for the magnetic field H and electric 
displacement field D such that 

/ —D x -DV —D z \ ( 



G 



a/3 _ 



D x -H z 
Dy H z 
y D z -H'y H x 



H'y 

-H x 



a/3 



—D x 

-D'y 

\ -D z 



D x 


Dy 


D z \ 





-H z 


Hy 


H z 





-H x 


-Hv 


H x 


o J 

(85) 
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Finally we introduce the polarization tensors which contain both the mag- 
netization M and polarization P 



yap 



( 



\ 






px 


py 


pz 


px 





-M z 


My 


py 


M z 





-M x 


pz 


-My 


M x 






\ 



) Pa/3 



( 





px 


—py 


-P z 




px 





—M z 


My 




py 


M z 





—M x 


\ 


pz 


-My 


M x 






(86) 

In terms of equations (83), (85) and (86) we can rewrite equation (2) and 
equation (3) as a single tensor equation 



Qa.fi _ pad _ ^pap 



(87) 



Furthermore, the equations appearing in (5) can be rewritten as a single 
tensor equation 



c • 



(88) 



In this the four current is defined as J a = (cp, J) and the subscript / refers 
to free currents and charges. In what follows we assume that the material 
under study does not contain any free currents or charges. The covariant 
Lagrangian density for the electromagnetic field can be written as 

Cfp = C-f + £f.j 

£ r = —p „F a/3 

16vr al3 

Cfj = — J Q A a 
c 



(89) 



In matter J a can be written as 



Ja = cd^Pt 



13a 



(90) 



This expression automatically satisfies the conservation law d a J a = 0, due 
to the antisymmetry of Pp a . The field equations (88) can be derived by 
taking the variational derivative of equation (89) with respect to the four 
potential A a . Writing C given in equation (89) in terms of the standard 
notation leads to the expression 

£ F p = —[E 2 - B 2 }+ P ■ E + M ■ B + V • ($P + M x A) + -8AA- P) (91) 
8ir c 

Notice the differences between equation (7) and the covariant result given 
in equation (91): 
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1. The covariant formalism contains the divergence term V- (&P+M x A) 
which is absent in the previous formalism. This term will contribute 
only on the boundary of the domain of integration and can be ignored 
if we assume that the domain of integration is over the entire space 
and the fields vanish at infinity. 

2. The covariant formalism contains the total time derivative -dt(A ■ 
P) which is absent in the previous formalism. This term will not 
contribute to the field equations and can be ignored. 

3. The coupling between the magnetization and magnetic field is differ- 
ent. In the covariant formalism the induction density B is coupled to 
the magnetization instead of the magnetic field H. 

A. 3 The energy momentum tensor of the electromagnetic 
field 

The energy momentum tensor in terms of arbitrary fields r\ p is given by 

T '- = m^ a ^- 5 " c (92) 

For the electromagnetic field, when we calculate the energy momentum ten- 
sor T@[FP] using Cfp, this results in 

T&FP] = ^-d a A^ + ApdaPrt - 5iC FP (93) 

The energy density is given by Tq[FP] which can be written as follows: 

T°[FP] = -!- [E 2 + H 2 — (4ttM) 2 ] + V • (-^ + A x M) (94) 

The divergence term V • (^ + A x M) will only contribute as a boundary 
term to the electromagnetic energy and can be ignored if we assume that 
the domain of integration is over the entire space and the fields vanish at 
infinity. Notice that for the covariant Lagrangian density the energy density 
is different from the result of equation (9): it contains a negative quadratic 
term in the magnetization instead of a positive one. If the magnetization is 
a given static quantity, this is just a constant that can be ignored. On the 
other hand, if the magnetization is a dynamic degree of freedom, this term 
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can render the energy negative. Of course, if the magnetization is dynamic, 
it will have its own Lagrangian which will have its own contribution to the 
energy momentum tensor; this will be discussed in the following subsections. 

(One would expect that the relativistic energy density in equation (93) 
would smoothly approach the non-relativistic expression in equation (9) as 
c — > oo. Apparently this is not the case, because of the different signs of 
the quadratic magnetization term. In truth, when one goes deeper, into the 
atomic theory of magnetization, one recognizes that both the orbital and 
spin contributions to the magnetization M are proportional to the Bohr 
magneton This vanishes in a theory where c — > oo.) 

Finally we derive the Poynting vector by calculating Tq [FP] . This results 

in 

T*[FP] = {±-ExH + h t {M x A - M) + y x (95) 

which is the "correct" form if we ignore boundary terms and total time 
derivatives. 
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